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Abstract 

Fibonacci chains are special diatomic, harmonic chains with uniform nearest neighbour in- 
teraction and two kinds of atoms (mass-ratio r) arranged according to the self-similar binary 
Fibonacci sequence ABAABAB A..., which is obtained by repeated substitution of A — ► AB 
and B — » A. 

The implications of the self-similarity of this sequence for the associated orthogonal polyno- 
mial systems which govern these Fibonacci chains with fixed mass-ratio r are studied. 


1 Introduction 


Fibonacci chains are linear diatomic chains with nearest neighbour harmonic interaction 
of uniform strength n and the two masses (ratio r = mi/m 0 ) follow the pattern of the binary 
sequence {h(ra)}i° obtained by repeated substitutions <r in the following way. 

<r(l) = 10 , <t(0) = 1 , (1) 

starting with 0. By definition cr(uv) = <r(u)<r(w) for any two strings u and v. <r n (0) = H n is a string 
of length \H n \ = F n+ i, where F n = F n . t + F n _ 2 , n = 2,3, .... F 0 = 0, F x = 1 are the Fibonacci 
numbers. h(n) is defined to be the n'th entry of the half-infinite string Hoc ;== lim n -*ooH n ■ E.g. 
Hs _ _ 10110101, h(l) - 1, h( 2) = 0, etc. (1) is called the Fibonacci substitution rule, and 

the masses of the half-infinite chain are taken to be 


This sequence {/i(n)}~ is self-similar because the string Hoc satisfies <r(Hoo) = H^ . Aperiodicity 
follows from this invariance, or fixed point, property. (This sequence is in fact also quasiperiod! c, 
but this does not concern us here.) 

Chains of this type have been considered as models of binary alloys [1]. For instance, one may 
consider chains with an elementary unit determined by the first N members of the (h(n)} sequence 
and repeat it periodically, using certain boundary conditions. This then corresponds to (AB) 00 
chains for N = 2, ( ABA)°° chains for N = 3, etc. 


The dual of such chains (with equal masses but two spring constants k 0 and «i following the 
pattern of the Fibonacci substitution sequence) are related to one-dimensional quasicrystals [2], 
One can also make contact to artificially manufactured superlattices [3]. 

Originally such Fibonacci chains were considered as models for the study of the regime in between 
periodic and random structures. [4, 5]. 

The purpose of this workjis to write down the identities which are satisfied by the characteristic 
polynomials of these Fibonacci chains due to the self-similarity of the substituioit sequence {h(n)} 
which determines the pattern of the masses of the oscillators. These identities will be expressed 
in terms of the 2x2 transfer matrices M n which are unimodular and real. The matrix elements 
are given by the characteristic polynomials (sl r )(x)}, where r is the mass-ratio of the two types 
of atoms and x is a normalized frequency squared (x - u) 2 /2u%, a % = /c/m 0 ). The zeros of 
■!>j^(x) determine the eigenfrequencies of finite Fibonacci chains with N atoms and both ends 
fixed. One also encounters so-called first associated polynomials {Sj^(x)}. They correspond to 
a right shift by one unit in the substitution sequence. Hence, the zeros of sj^(x) produce the 
eigenfrequencies of chains with masses rn^ 2 ) = m-o, ■••y‘ rr ^h(N+i)- Both r— families of polynomials 
generalize Chebyshev’s {S n (y)} polynomials (5_ x = 0, 5 0 = 1, 5 n = yS n _!-S n _ 2 ) to two variables 
with the identification 

SW*) = $?>(*) = S„( 2(1 - *)) . (3) 

They constitute, for fixed mass-ratio r, systems of orthogonal polynomials and have been studied 
in some detail in refs. [6, 7, 8, 9]. 

2 Fibonacci Chain Polynomials 

For the Fibonacci chains («,m^( n )) defined in section 1 the equation of motion for longitu- 
dinal, time-stationary vibrations q n (t) = q n exp (icot) are 

?n+ 1 + q n - 1 - Y(n)q n = 0, n = l,2,... ( 4 ) 

with 

Y{n) = 2(1 - u> 2 /{2u> 2 n )) , = «/m Mn) . (5) 

We use the two variables r = mi/m 0 and x = u> 2 /(2u>Z). We put Y{n) = Y if h(n) = 1 and 
Y(n) = y if h(n) = 0. Hence 


Y = 2(1 -rx) , y = 2(1 - x) 


( 6 ) 


The equations of motion are rewritten with the help of the SL(2,R) transfer matrix R n : 


(V)-‘(;.)-(T 



( 7 ) 


Rn is either R\ or Rq depending on the F^n) value, i.e R„ = For the half-sided infinite 

chains considered here iteration leads to 



( 8 ) 
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with the inputs qi and q 0 (the mass at site number 0 is irrelevant). M n is real and unimodular. 
The recursion M n = RrMn - 1 with input Mi = Ri leads to 


M n = 



-s n - 2 ) ' 


( 9 ) 


where the recursion formulae for the generalized two-variable Chebyshev polynomials are 


S n = y(n)5 n _ x - S n _ 2 ^ , 5_! = 0, S 0 =l , (10) 

5„ = Y(n + l)5 n _i — S n -2 , S- 1=0, 50 = 1 . (H) 


These polynomials generate certain combinatorial numbers [10]. The meaning of these numbers 
can be understood if one uses the intimate connection of the Fibonacci substitution sequence with 
Wythoff’s A and B sequences 


n — 1 

A(n) = n -f ^^h(fc) , B(n) = n + A(n) . (12) 

fc=i 


These sequences {A(n)}?° and {£(n)> cover the positive integers in a complementary way: every 
number N > 0 is either an A - or a B-number. For an A-number n (i.e. n = A(m) for some 
m) h(n) = 1, and for a B — number n (i.e n = B(m) for some m) h(n ) = 0. Wythoff s sequences 
are a special case of Beatty sequences: A(n) = \nip J, B(n) = L^^J) will 1 = + 

golden mean. 

The characteristic polynomials {5i r ^(x)}, obtained from {5 n (Y, y)} by replacement of 1 and 
y according to eq.(6), constitute, for fixed mass-ratio r, a system of orthogonal polynomials. 
(5i r l(a:)} are the first-associated orthogonal polynomials. 


3 Self- Similarity Identities 

The string, or ’word’, H <*, defined in section 1 is invariant under the inverse substitution a , 
with <r -1 (l) = 0, cr _1 (10) = 1. This is equivalent to the self-similarity of the sequence {Mu)}i° 
which is shown in the FIG. 



FIG. Self-similarity of the sequence {^(w)}J° . Circles stand 
for th? value 1 (A-numbers n), and disks stand for the value 
0 (J3-numbers n). (T _1 (10) = 1 , cr -1 (l) = 0 . Level (/) is 
mapped to level (/ + 1) by <r _1 . 
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The upper level, called (/) in the FIG., shows the numbers marked as A— and 0— numbers. The 
h(n) value is 1 or 0, denoted by a circle or disk, respectively. When the substitution <r -1 is 
applied one reaches the next higher level, called (l + 1) in the FIG., on which the same sequence 
is reproduced. Let the position of the n' s number at level (/) be for 1 = 0, 1,.... Level l = 0 is 
assumed to correspond to the original sequence. Then one finds for p £ No = N U{0} 


x 

x 


0 + 1 ) 

Mp) 

O+i) 

b(p) 


.<0 

X B(p) i 

JO 

X AB(p) 


(13) 

(14) 


where AB(p) stands for the composition A(B(p)) of Wythoff’s sequences. E.g. The number 
.4(4) = 6 at level (l + 1) occurs in the FIG. at the same position as B{ 4) = 10 at level (/), or 
.0(2) = 5 at level (l + 1) corresponds to .40(2) = .4(5) = 8 at level (I). 

Iteration, depending on the parity of the level number, leads to 


(2fc+l) _ (0) 

X A(p) ~ X B*+‘(p) 

} 

J 2 *) _ (0) 

X Mp ) ~ X AB“(p) » 

(15) 

(2*+l) _ (0) 

X B(p) ~ X XB‘+*(p) 

} 

J 2 *) - J°) 

X B(p) - X B k + 1 (p) ■ 

(16) 

Consider the level (f + 1) transfer matrix 




M< i+1) = 

h(n) 

" mR h{ 1) » 

(17) 

satisfying the recursion relation 




<■> = , 

^(°) 

n 0 

■"■■(! V) ■ 

(18) 

RH* l) = R? , 

r[ 0) 

■ - - fi V) - 

(19) 

Iteration leads , with M n = Mn \ to 




at" = m F( w> 

, at" = , 

(20) 

with the Fibonacci numbers F n . 
Due to (15) and (16) one has 




¥am 1) = 

> 

M A(p) ~ M AB k (p) , 

(21) 

KT = 

1 

*4(3 = m b ^(p) • 

(22) 

The recursion at each level is 




Mt l) = aQ$ , 

mJ ,+ 1) = 0< i+1) . 

(23) 


Combining iteration and recursion, in a systematic way, leads to transfer matrix identities for 
level (0), i.e. for the original matrices M n of eq. (8). One finds alltogether six families of such 
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identities, depending on the parity of the level one starts with and the specification of the index. 
These identities are, for m G N and k G N, 


(I) 


(II*) 

M B k A(A(m)+\) 

(Hb) 

M B k A ( B ( m ) + l) 

(III) 


(IV a) 

M AB k A(A(m)+\) 

(IVb) 

M AB k A{B(m)+l) 


= AfFjt+i ^AB k A(m) ) 
= -^B*+ l (m) > 

= Mf 3Ic+1 M AB k-n( m ) 

= Mgk A ( m ) , 

= M Fl(k + l) M AB k + '{m) 

= ^(t + i) ^B*+ 3 (m) 


(24) 


(/), (//o), (776) and (777) result from odd levels / = 2k + 1, with n put A7?(m), AA(A(m) + 1), 
AA(B(m) + 1 and £(m), respectively. (777), (7Va) and (IVb) result from even levels l = 2fc, 
with n put B(m), A(A(m) + 1) and A(B(m) + 1), respectively. 

E.g. (7) and (777) produce for m = 1, due to >1(1) = 1, £ fc+1 (l) = F 2fc+3 and AB*(1) = T^+j), 
identities which are the well-known recursion formula for transfer matrices with neighbouring 
Fibonacci number indices 

M Fn+1 = M f „_, M Fn . (25) 

Not all eqs.(24) are independent. E.g. if one puts m = B(p) + 1 in (7), replaces k by A: + 1 and 
combines it with eqs. (IVb), with A: — » k — 1 and m — > p, one finds eqs. (77 a), due to the identiy 
£(p) + 1 = i4(j4(p) -1- 1) and eq.(25) for even n. However, eqs. (77a) provide identities for M B k A{p ) 
which complement those obtained from eqs. (776). 

It is possible to combine (7) of eq.(24) with (777) specialized to m > j4(jti) and use (7) again 
with k k — 1 and m — ► A 2 (m). Continuing this process one finds for A: G N and m G N 


(7') = M Fjk+l M Flk • • ■ Mf 2 M BA ik( m ) 

(in') M AB k (m) = Mp 3k M -M Fi M BA ik- i( m) (26) 

(7) and (III) in (24) can be replaced by both eqs. (26), and the other eqs. of (24) can be rewritten 
using (26). 

The sum of the indices of the transfer matrices on the r.h.s. of eqs. (24) and (26) have to match 
the index of the l.h.s. This fact produces families of identities among iterated Wythoff A and B 
sequences. A detailed investigation of these Wythoff composites identities will be given elsewhere. 
All of these identities can be rederived as corollaries of a new theorem relating two seemingly dif- 
ferent unique number systems: the Wythoff- and the Zeckendorf- (or Fibonacci-) representations. 

The transfer matrix identities (24) are equivalent to those for their matrix elements, i.e. the 
characteristic polynomials {5 n (h r , t/)} and {^(IK, j/)}. In order to derive them one rewrites the 
indices of all matrix elements as Wythoff composites. Consider, for example, (7). For the elements 
of one employs the simple identities B k+1 (m) — 1 = B(B k (m)) — 1 = A 7 B k (m) and 

7? fc+1 (m) - 2 - ABAB k ~ 7 (m). The last identity can be proved for m = A(p) and m = B(p) 
separately. On the r.h.s. of (7) one rewrites the indices of the matrix elements with the help of 
the identities T^k+i = B k (l), F 2 *+ 1 — 1 = A 2 B k ~ 1 ( 1), F 2 k+i — 2 = ABAB k 2 (1) for k = 2,3,..., 
and 7^3 — 2 = 0. Moreover, AB k A(m) — 1 = BAB k ~ L A(m), AB k A(m) — 2 = A 3 B k x A(m). 
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Finally, (/) decomposes into the following four sets of eqs. 


(MM)) 

(MM)) S A 2 B k( m ) 

(M 2 >1)) S A 2 B k (m) 

(M 2 > 2 )) S A BAB k ~ i (m ) 


$B k ( 1) ^AB*j4(m) — *^BXB*-M(m) ) 

*^B*( 1) S. BAB k - i A{m ) ~ ^i J B*-»( 1) ^A 2 B k ~ l A(m) » 
S.4 J B*- l (l) — ■S’jtBXB*-^!) S fl XB*- l >l(m) > 

^ A 2 B k ~ l (l) ‘^BAB k ~ 1 A(m) ~ $ ABAB k ~ 2 (l ) S A »B k - 2 A(m) 


The last two sets of eqs. hold only for fc = 2,3, .... For k = 1 one has 


(27) 


S A 2 B(m) — Y SABA(m) ~~ , 

S ABA(m ) = Y *^A 4 (m) • (28) 

The other eqs. in (24) decompose in a similar way. The arguments of the polynomials is always 
(F, y) , which can be replaced using eq.(6). 

This concludes the derivation of the self- similarity eqs. for the Fibonacci chain polynomials. It is 
clear that further work is needed in order to extract from this gamut of eqs. information pertaining 
to chain properties, like structure of spectra and displacements. 
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